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Abstract 

In order to explain neutrino oscillations, it is usually assumed that neutrinos have tiny 
masses, and that the mass eigenstates are superpositions of flavor eigenstates. However, the 
nature of neutrino mass and the mechanism by which flavors mix have not been well 
explained. In this article, we show that the nature of neutrino mass is a potential energy 
inversely proportional to the neutrino energy, which origins from the exchange of 
experimental neutrinos with cosmic neutrino background. The essence of neutrino flavor 
transition in neutrino oscillations is the exchange between neutrinos of different flavors. 
With the assistance of data obtained from neutrino oscillation experiments, the potential 
energy matrix is determined to contain only one independent parameter. The elements of 


the potential energy matrix have characteristics consistent with the distribution of 


background neutrinos predicted by the standard cosmological model. Our model can not 


only explain the normal neutrino oscillations, but also provide a unified perspective to 
explain the long-standing anomalies in some neutrino oscillation experiments. Further 


experimental study of these anomalies will provide a feasible way to test our model. 


Introduction 

It has been experimentally demonstrated more than twenty years that neutrinos can change their 
flavor periodically with distance over energy while they propagate, which is known as neutrino 
oscillations!?. In order to explain this amazing phenomenon, it is usually assumed that neutrinos 
have tiny masses, and that the mass eigenstates are superpositions of flavor eigenstates®’. In 
quantum mechanics, the waves associated with neutrinos of different mass have different 
frequencies. This gives rise to time-dependent quantum-mechanical interference effects, and the 
composition ofthe beam will change along its path. The existence of neutrino oscillations challenges 
two assumptions in the Standard Model of Particle Physics (SM): neutrinos have zero rest mass and 
the flavor of leptons is conserved??. In general, the mechanism of neutrino mass would require the 
addition of particle content to the SM that has never been observed, while the interpretation of 
leptonic flavor mixing usually requires the introduction of special flavor symmetry groups!?-?. 

Despite strenuous efforts, a comprehensive understanding of the origin of neutrino mass and 
flavor mixing remains elusive. It is worth noting that, there is still no direct evidence to date for 
neutrino mass!3, and no other physical process violating lepton flavor conservation has been 
identified!^!6, Furthermore, considering the long-standing anomalies observed in certain neutrino 


oscillation experiments that still lack satisfactory explanations!"?!, it seems that we need to revisit 


our understanding of neutrino oscillations. In this regard, a theoretical model in condensed matter 
physics can give us some inspiration. In some magnetic materials, a moving electron can exchange 
with other electrons it encounters, which is used to explain the ferromagnetism exhibited by these 
materials (Itinerant Electron Model)?. Although Coulomb forces do not exist between neutrinos as 
they do between electrons, the exchange between neutrinos itself can explain the neutrino 
oscillations, as I will show below. 

The explanation is based on a fundamental quantum mechanical principle - the exchange effect 
between identical particles. It should be emphasized that the exchange effect described here is not 
the neutrino-neutrino scattering mediated by the electroweak interactions. Even in the absence of 
interaction transmitted by vector bosons, neutrinos still adhere to certain statistical laws (neutrinos 
are fermions and follow Fermi-Dirac statistics). This fact indicates that when the wave functions of 
two neutrinos overlap in space, there must be an exchange effect between them. The situation here 
is analogous to the exchange of electrons between two hydrogen atoms in close proximity to each 
other when forming a hydrogen molecule. The Standard Cosmological Model (SCM) predicts that, 
the Cosmic Neutrino Background (CvB), a relic from the early university when it was about one 
second old, consists today of about 112 cm? neutrinos plus antineutrinos per flavor?*?, It is the 
largest neutrino density at Earth?°. If neutrinos were massless, the CvB would be blackbody radiation 
at T= 1.95 K = 0.168 meV, corresponding to a de Broglie wavelength of about 6 millimeters. Any 
neutrino detected in the neutrino oscillation experiments, which we might call experimental neutrino, 
actually travels through the background neutrino sea before being detected. The exchange of 
experimental neutrinos with background neutrinos is the essence of neutrino oscillations (Fig. 1). 
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Figure 1. Schematic of neutrino oscillations. When an experimental neutrino (the red one) travels 
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through the background neutrino sea, it can exchange with the background neutrinos (the green 


ones), which leads to the neutrino oscillations. 


Exchange effect and potential energy matrix 

Consider an experimental neutrino ve of energy Ei exchanging with a background neutrino ve of 
energy E2. We denote the background neutrino by v, (a = e, u, x), which is used only for the sake of 
narrative convenience, and does not imply any essential difference between the background 
neutrinos and the usual experimental neutrinos. For the experimental neutrino ve, it can be in either 
the Ei or E» state due to the exchanged with the background neutrino ve. This can be viewed as a 


two-state system. Assuming that the exchange amplitude is A, the corresponding Hamiltonian is 


(E A 
H=| 2 (1) 
A E, 
Its two eigenvalues can be derived as (See Part A in SI for details) 


E' = (E, + E; T4 (E1 — E2)? + 44? ) /2, (2a) 


E',= (Œ +E, —Jf (Ey — Ey)? + A4?) /2. (2b) 

Assuming (Ei- £2)? > A? and Ei > E», we have 
B= E ATE (3a) 
E', = E; — A! JE, (3b) 


According to Eq. (3a), the neutrino ve with energy Ei acquires a potential energy inversely 
proportional to its energy by exchange with the background neutrino ve. In explaining the neutrino 
oscillations, it is usually assumed that the neutrino has a tiny mass m, so that the time-dependent 
phase factor of the wave function can be approximated to (in natural units A ^ c = 1) 

exp(iEt) = exp(i./ p? + m2t) ~ exp[i(p + m?/(2p))t]. (4) 
where p is the magnitude of the neutrino momentum, and p > m is assumed. For the neutrino with 
zero rest mass, Eq. (3a) can be rewritten as Ej ~ p + A*/p. Compared to the phase factor in Eq. 
(4), it can be seen that in terms of the effect on the phase factor of the wave function (which is the 
key to explaining the neutrino oscillations), the fact that the neutrino has a tiny potential energy 
inversely proportional to its energy (V = A’/E)) is equivalent to giving the neutrino a tiny apparent 
mass m, as long as m = V2|A| is taken. 

As a neutrino ve travels through space, it can exchange with more than one ve at the same time. 
We write down the total amplitude as Ace and the corresponding potential energy as Vee. We further 
assume that the neutrino ve can also exchange with v. and v; in the CvB. The flavor of the leptons 
may be a property similar to spin polarization, and just as photons with different polarizations are 
still identical particles, neutrinos with different flavors can be identical particles in some sense, and 
thus have exchange effects as well. In contrast, charged leptons of different flavors are not identical 
particles because they have different rest masses. The exchange of ve with v. orv- has two effects: 
one is to give the neutrino the corresponding potential energy Ve. or Ver, and the other is to change 
the flavor of the neutrino. Once the flavor is changed, e.g., when ve change into v, (v;), the v, (v2) 
can also be exchanged with ve, vy, and v: to obtain the corresponding potential energies Vie (Vre), 
Vim (Vru), and Vir (Ver). Phenomenologically, the neutrino va transitions between different flavors as 
it propagates and has a potential energy corresponding to each transition process. This makes the 
potential eigenstates of experimental neutrinos different from the flavor eigenstates. Considering 
that for flavor conversion, the exchange between va and vg (a, D = e, u, v; a Z B) and the exchange 
between vg and v, are two opposite processes, i.e., the former transforms v, into vg while the latter 
transforms vp into va in terms of the outcome, we can write the neutrino's potential energy matrix in 


flavor picture as follow 


Vee Veu — Vie Ver — Vre Aĉe Acu z Aiie Acr A Aie 
9 =| Veu — Vue Vu Yur — Veu | = : Abu — Aiie Ain Ajit — Atu » (5) 
Ver — Vre Var Vru Ver Acr — Ate WAS a. 


where £ is the energy of the experimental neutrino and Aag is the exchange amplitude corresponding 
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to Vag. According to the SCM, the three background neutrinos have similar distributions, with v, 
and v: having nearly the same distribution and ve having a slightly different distribution???^. We 
expect such a characterization to be reflected in the matrix elements of V. In addition, we expect 
Aa and Aba (a Z D) to be slightly different. The exchange of vg and va (a + B) would add an extra 
vp to the background neutrinos, and the process is subject to the restriction of the Pauli exclusion 
principle, whereas there is no such restriction when va and v, are exchanged. 

From the above analysis, the eigenvalues of the matrix V determine the three apparent masses 
of the neutrinos, while its eigenstates determine how the neutrino flavor is mixed. Define the matrix 


A’ as 


Ace Aen Aer 
Â? =| Aje Auu Am), (6) 
Ate Au Ate 


where, Aba = Aga (& = e, u, t); Abu = Age = Ay — Age» Ap = Atu = Age — Az,» and 


A2, = A2, = AZ, — Aze. From Eqs. (5) and (6), we have V = A?/E, and thus there is a relation 
V; = A?/E between the eigenvalues Vi of V andthe eigenvalues A? of A?. The three eigenvalues 
of the matrix A? can be solved as 


A? = 4? + 2|S|!/3cos[(@ — 41)/3], (7a) 
A2 = 4? + 2|S|'/3cos[(0 — 21)/3], (7b) 
A2 = 4? + 2/5|!/3cos(0/3); (7c) 


where, |S| and 0 (0 x 0 <7) are functions of the matrix elements of matrix A?, and A? = (A2, + 
Az, + 44)/3 (See Part B in SI for details). 

According to the previous analysis, the three eigenvalues A? (i= 1, 2, 3) correspond to the three 
apparent masses Mi (i = 1, 2, 3) of neutrinos, and satisfy M? = 2A?. From Eqs. (7a-c) we have 


M2 — M2 = 2(A3 — A?) = 4N3|S|!"sin(21/3 — 0/3)) > 0, (8a) 
M2 — M3 = 2(42 — 42) = 4N3|S|!"sin(z/3 — 0/3) > 0, (8b) 
M2 — M2 = 2(A3 — A2) = 4N3|S|!/?sin(0/3) > 0. (8c) 


Current neutrino oscillation experiments have not been able to fully determine the ordering of the 
three neutrino masses mi (i = 1, 2, 3). There are two possibilities, which are referred to as normal 
ordering (NO, m3 > m» > mı) and inverse ordering (IO, m» > mı m3), respectively. From Eqs. (8a- 
c), we have mı = Mi, m= M», ms = Ms for NO, and mi = Mo, m» = Ms, ms = M; for IO. Either 
in the case of NO or IO, using the available data from neutrino oscillation experiments, we can 
determine the values of |S|'/? and 0. Taking the case of NO as an example, using Eqs. (8a) and 
(8c), and taking the best-fit values" 42, = m2 — m? = 7.42 x 10 ^eV?and 42, = m2 — m? = 
2.514 x 10 ?eV?, we obtain 

sin(0/3)/sin[(2z — 9)/3] = (m2 — m?)/(m$ — m?) « 0.0295. (9) 
Solving the above equation gives 0 = 4.46°. Substituting it into Eq. (8a) gives |S|1/3 ~ 
4.130 x 107* eV?. From Eqs. (7a-c), we have 


A? = 4? — 4.314 x 107 ^eV?, (10a) 
A2 = 4? — 3.943 x 10-*eV?, (10b) 
A2 = 4? + 8.257 x 107^eV?. (10c) 


To explain the neutrino oscillations, it is usually assumed that the flavor eigenstates |v4) (a = 
€, I, t) of neutrinos are a superposition of three mass eigenstates |v;) (i = 1,2, 3), i.e., 
Iva) = Xi-i Uailvi). (11) 
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The U4; are elements of the unitary matrix Ü, which is known as the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) lepton mixing matrix to be parameterized in terms of three mixing angles (012, 423, 


013), and one phase parameter dcp, 


—iô 
C12€13 5412013 S430 "CP 
jj — id id, 
U = | —5412€23 — €12513323€ ^^" — €12€23 — $12543523€'°? €13823 b (12) 
ió id, 
512523 — €123413C23€ €P  —0415323 — $412513C23€0 "^" C13C23 


where, sj— sin, cij = cosOj. From the above analysis, m? = 242 (in the case of NO), so the 
eigenstate of the apparent mass mj is the eigenstate of A2. This means that in the flavor picture, the 
eigenstate corresponding to the eigenvalue A? is the i-th column of the matrix Ü. We can 
diagonalize the matrix A? (Eq. (6)) with matrix JU, that is, 

Ut AU = diag (4?), (13) 
where the diagonal matrix element in row i of the diagonal matrix diag(A”) is A2. The Eq. (13) can 
be equivalently written as 


A? = Udiag (4?)0*. (14) 
From Eqs. (12) and (14), the elements of matrix A? can be solved as 
A2, = A? ch ci AA + $2,c?,4A3 + 57,AA3, (15a) 


A2, = (761231201323 — C150135135823€ CP) AAT 
t (c12512€13€23 E 5220435138236 CP )AA? (15b) 
0435138256 CP AAS, 
A2. = (€12512€13523 — Cf2C13S13C23€ CP AA? 
t (—612512€13523 = S25C13513C23e 9c" JAAS (15c) 
0435430236. CP AAZ, 


2 _ 2 id 2 
Ane = (—C12812€13C23 — €12013513823€ ° )AA1 


t (c12512€13€23 = 5220435138236 CP )AA2 (15d) 
04351382369? AAZ, 
Aju = A2 + (s2,c33 + c] 523533 + 2012512513C23523C0SOcp) AA? 
+(cf2C33 + SfoSf3833 — 2€42812513C23823COSScp AAS (15e) 
tchs2,AA2, 
Af. = (—S$2C23823 + C12513023823 + €42512513C33@ ‘OCP — c455,55,352,0 CP AAT 
t (—c25623523 + S25525€23523 — Ci2512513C246 cP + 0125425435236 c )AA2 (15f) 
-C13C235234A$, 
A2, = (€12512€13523 — C12€13513 023€ OC? )AA1 
t (—612512€13523 = S2, C13543 C236 cP )A42 (15g) 
0435430236 P AA, 
Az, = (—512C23523 + Ci251323523 + C42512513 0256/9? — c4551554352, e OCP AAP 
t (—c35623523 + S25525C23523 — Ci25125413 C236 OC? + 125425438236 (Oc? AAR (15h) 
+C73C235234A3, 
A2, = A? + (sis, + c]5sisc]s — 2¢12512513C23523C0SScp) AA? 
+(cf2833 T S12513C23 + 2€1251251323523C0SÓcp)AA7 (151) 
+e 030A, 


where AA? = 2|S|*/3cos((@ — 41)/3) , 443 -—2]|S|!?cos((0 —21)/3) , and 443 = 
2|5|!/3cos(0/3). 

In the case of NO, the best-fit value of dcp is 195°, and here we may wish to take dcp = 180? to 
simplify the analysis and to be consistent with the earlier taking of the matrix A? (Eq. (6)) to bea 
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real matrix. In principle, for any dcp, we can calculate the matrix A? with the help of Eqs. (15a-i). 
However, if the matrix elements of Ü are imaginary, then the non-diagonal matrix elements of A? 
will also be imaginary. This means that the matrix V (Eq. (5)) needs to be introduced with the 
corresponding phase factors in its non-diagonal matrix elements. For the three mixing angles, we 
take the best-fit values?’ 412 = 33.44°, 023 = 49.0? and 013 = 8.57". Substituting the above taken values 
and Eqs. (10a-c) into Eqs. (15a-i), with 1074 eV? as the unit of the matrix elements, we have 


_ (44-3925. -1275 1.332 
A =| -1275 42-2825 5956 | (16) 
—1.332 5.956 — 474 1.100 


The parameters (A?;(= m? — m?), Oj, Ocp) measured by the neutrino oscillation experiments 
have completely determined the values of Aab (a + B) and A24, — Age (Eqs. (15a-i)). This can 
also be seen directly from Eq. (16). Also, it can be seen that Aab (a + B) and A2, — AR g are in 
the same order of magnitude. This is consistent with our expectations. It can also be seen that 
between the diagonal elements A2,, A7 
closer. All these features are still present when taking dcp = 0 and even in the case of IO (dcr = 0, x) 


(Table 1, see Part C in SI for details). 


u and AZ, the values of Aj, and A2, are somewhat 


Mass Order | àcp A2, Abu AA A2, A2 A2; 
NO m | À4-3.95| A?+2.825| A? «1.110. -1.275 | -1.332 | 5.956 
NO 0 | 4?-3.95| A?+2.775| A? «1.496| 1.496 | 1.077 | 5.963 
IO m | A?+3.868| A? —2.841| A? -1.028|. 1.524 | 1.088 | -6.164 
IO 0 | A? +3.868| A? -2.891| A? —0.977| -1.304 | -1.344 | -6.157 


Table 1. The calculated values of matrix elements of A? (in the unit of 10^eV? ). 


The experimental parameters Aj and the mixing matrix Ü (in terms of 6; j) of neutrino 
oscillations can be derived from the matrix A? (Eq. (16)), which suggests that these parameters are 
in fact determined by Aab (a + B) and A24, — Apps i.e., that the neutrino oscillations are mainly 
determined by the absolute differences between the three background neutrino distributions. The 
relative differences between 42, (a =e, u, T) (ie., |A24 — A$gl / A2 ) depend on the value of 
A2. If A? — 500 is taken, the relative difference between A2,(a = e, u, T) ison the order of one 
percent by Eq. (16); if A? — 5000 is taken, the difference is on the order of one thousandth. From 
the above analysis, the study of neutrino oscillations provides information on the differences in the 
distribution of background neutrinos with different flavors. 

Understanding the anomalies in neutrino oscillations 

Given A? in Eq. (16), the matrix A? is uniquely determined. Although the available 
experimental data on neutrino oscillations does not determine the value of A2, the anomalies in 
some neutrino oscillation experiments provide some hints as to the range of possible values of A42. 
There are four long-standing anomalies in the short-baseline neutrino experiments! 7?!282?, Two 
arise from the apparent oscillatory appearance of electron (anti)neutrinos in relatively pure muon- 
(anti)neutrino beams. The other two anomalies are associated with an overall normalization 
discrepancy of electron antineutrinos expected both from conventional fission reactors (the Reactor 
Antineutrino Anomaly) and in the radioactive decay of ’!Ga (the Gallium Anomaly). In these latter 
cases, no oscillatory signature is observed, but the overall normalization deficit can be ascribed to 
rapid oscillations at a high Aj; ~ 1eV? that are overaged out and appear as an overall deficit. The 


neutrino experiments in which the above anomalies occur share several common features: firstly, 


6 


they all involve high-flux neutrino sources. For example, in the case of the Ga anomaly [19], the 
neutrino flux ®, near the neutrino source is on the order of 101* cm~*s7}; secondly, there are 
obvious tensions in the range of possible values of A (and 6;;) obtained from different 
experiments; thirdly, a clear oscillatory signal dependent on L/E (L is the detector-to-neutrino-source 
distance, and Æ is the neutrino energy) has not been observed in almost all the cases (except for the 
Neutrino-4 experiment?) In the following we will show that our model provides a unified 
perspective for understanding the neutrino anomalies mentioned above. 

In the previous analysis, we implicitly assumed that the CvB is not perturbed by exchanges with 
experimental neutrinos, and thus the matrix Â? is definite. This may not be true for the short- 
baseline neutrino experiments described above. The SCM predicts that the flux of the cosmic 
neutrino background is ®, ~ 101? cm~*s71. If the neutrino flux ®, near the neutrino source used 
in a neutrino oscillation experiment satisfies ®, >> ®,, then the CvB near the neutrino source may 
be disturbed by exchanges with experimental neutrinos. This will change the values of Aab (a + 
B) and Aža — A$g- As a result, the mixing matrix Ü and the apparent mass m; of the neutrinos 
are changed. At a given perturbation, Aj; may increase, and this will lead to more rapid oscillations. 
The result would be a further enhancement of the perturbation to the background neutrinos. This 
nonlinear effect may allow Aj; to increase rapidly to the order of A?. We can make a rough 
estimate. In unit of 10 ^ eV?, taking A2, = 6000, AZ, — 5000 and A2, = 4000 (i.e., a 20% 
perturbation with respect to the mean value A? = 5000, specifically, 20% of v, transitions to v2) 
and simply taking A25 = Aga — Afp» ie. AZ, = 1000, A7, = 1000, and A2, = 2000. Direct 
calculations give the three eigenvalues of A2 ~ 7895, A2 ~ 4397, and A7 ~ 2708, respectively. 
(It is worth pointing out that, for experimental neutrinos with energy E of about 1 MeV? or greater, 
the potential energy (V; = A? / E) corresponding to the above values is quite small.) We have A2,— 
2(A3 — AZ) 10357. This means that if A? = 5000, a 20% perturbation could lead to 
hv 1 eV?. Because the details of the neutrino sources are different in various high-flux neutrino 
experiments, the perturbation of the CvB is also different. This explains why there is a tension 
between the fitted data from different experiments. In addition, the neutrino flux varies with distance, 
and thus the perturbation of the CvB also varies with distance, which causes Aj to vary with 
distance. It is anticipated that the oscillatory signal dependent on L/E will be eliminated, so that the 
experimental results reflect only some kind of averaging effect. From the above analysis, the study 
of anomalies in neutrino oscillation experiments will provide information on the absolute magnitude 
of the background neutrino flux. Such studies also provide a feasible way to test our theory. If the 
correlation between neutrino flux and the anomalies is systematically studied by experiments and a 
clear relationship between them is confirmed, this will be a strong support for our theory. 

The SCM predicts that neutrinos and antineutrinos in the cosmic neutrino background have the 
same distribution, which means that the results of the above analysis apply to both neutrinos and 
antineutrinos. On the other hand, if one analyzes the oscillation data for neutrinos and antineutrinos 
separately and confirms that Aj and 0j; have different best-fit values in the two cases, then 
according to our model, this would suggest that neutrinos and antineutrinos in the CvB have different 
distributions. As a final note, since the zero-mass neutrinos and antineutrinos possess opposite 
helicity, an exchange between them would violate the conservation of angular momentum and thus 
need not be considered. 

Conclusions 


In summary, we put forward a new interpretation for the neutrino oscillations. According to our 
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interpretation, the nature of the neutrino mass is a potential energy inversely proportional to the 
energy of the neutrino, which is acquired by exchange with CvB. The observed change in neutrino 
flavors in oscillation experiments results from the exchange between experimental neutrinos and 
background neutrinos with different flavors. We have analyzed the potential energy matrix of the 
experimental neutrinos using the parameters measured by the neutrino oscillation experiments, and 
the obtained results have features consistent with the distribution of the CvB predicted by the SCM. 
Moreover, according to our interpretation, if the neutrino experiments involve a neutrino flux much 
larger than that of CvB, the neutrino exchange will likely perturb the distribution of CvB, causing 
the apparent mass of the experimental neutrino to change. This provides a unified new perspective 
for explaining the anomalies in various neutrino experiments which usually involve high-flux 
neutrino sources. 

Unlike many existing neutrino theories that are difficult to test experimentally (e.g., the Seesaw 
Model?!), our model can be confirmed or falsified by feasible experiments. For instance, studying 
the relationship between the anomalies mentioned above and the magnitude of the neutrino flux 
would provide a feasible way to verify our model. Furthermore, our model suggests that neutrino 
oscillations do not imply that neutrinos have mass or that the conservation of leptonic flavor is 
violated. If future cosmological observations? show that the upper limit on the sum of neutrino 
masses (m; + m» + m3) is less than 50 meV, this would disprove the popular interpretation of neutrino 
oscillations, and thus be a piece of evidence in favor of our model. On the other hand, if experiments 
confirm the existence of neutrinoless BB decays or other lepton flavor violating decays, that would 
be a piece of evidence against our model. 

From the theoretical point of view, it would also be a strong support for our model if one could 
calculate the matrix (Eq. (16)) starting from the basic principles of quantum mechanics and the 
distribution function of background neutrinos predicted by the SCM. Further research is needed on 
related issues. In addition to the cosmic neutrino background, the entire universe is permeated by 
the cosmic microwave background. Just as experimental neutrinos are exchanged with background 
neutrinos, photons from distant galaxies are expected to be exchanged with the cosmic microwave 
background as it travels through cosmic space. The effect of this exchange is an interesting subject 
for study. 


Data availability 
The datasets used and analysed during the current study available from the corresponding author on 
reasonable request. 
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Part A. Find the eigenvalues of matrix H (Eq. (1) in the main text) : 


a [E A 
H= , 
A IE, 
Solution: 
Denoting the eigenvalues as. E/(i 21,2) , we have 
E, — E! A 
det ,|=0. (A1) 
A B-E 
That is 
(EY - (E, + E)E; + EE, - A =0. (A2) 
The solutions are 
El, =[E, + E, (E, - E,} +44°]/2. (A3) 


If (E-E, >A’ and E, E,, we have 


Ej, = [E + E, £(E, - EJ0* 24 / (E, - E»))/2 (A4) 

namely, 
ERE +A /(E-E,) =x E +A TE, (A5) 
E, x E, -A I(E- E,) = E,- A /E,. (A6) 
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Part B. Find the eigenvalues of matrix Ae 


(Eq. (6) in the main text) : 


2 2 
A; 77 A; 

2 2 
Ay, A, T f? 

2 2 
a A; 


where, A2, = Az, Aj, =A, and Aj =A’. 


ue’ uT 


Solution: 


Denoting the eigenvalues as A? (i 21,2,3), we have 


2 2 2 
A; a A; A 


2 2 2 

det] A, — AL-A4 
2 2 
Ay, Ai, 


That is 
(A?) — (AZ + a. + ANAY 


HALGA AU ACA SUP SA 


ee^ uu y TT ee^ TT e 


2 242 2 242 2 2: <2 2 2 2 2 2 2: 
FA (A) T Azi (Aa) st A. (A) E ACA A. E 2A, Arc Aer E 0. 
Let 
a--(A) + Aun +A’), 
2 2 2 2 2 2 2 x2 2 42 2 x2 
b- ALAL +A. Ar +A Ap (AY - (AL - (AZ, 
2 2 42 2 2*2. 2 2 2 2:42 2 2 2 2 
c= A (Aj) F A (A) + A. (A,,) E AA AS F 2A A A. £ 
We have 


2 
A; T 


y= (ARYA? 


u uT 


(A?) +a(A4 Y +bA? +c=0. 


The solutions are 


A =A? eod-gI2 de IA P121 +0? d-gl2- fg /44 f° /27 


A =A? o cerae g a PII S od|-el2- Jah 1a P127 


A =A +il-g (2492/44 Pim esa giae fm 


i2z/3 
, 


where, w=e and 


A! 2-al3- (A + A2, +A2)/3, 
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(Bl) 


(B2) 


(B3) 
(B4) 


(B5) 


(B6) 


(B7) 
(B8) 


(B9) 


(B10) 


f =b=a' 13 
= [AZ Ay, +A? A. + A> A> (ALY -(ALY -(ALy]/3 (B11) 
2 42 242 2.42 
SUAE EA OQ (43) 
g =c—ab/3+2a?/27 
2 2 2 2 22 2 2-82 2 D UND 
=-2A> A7 A? + 2A} (Az) * ALCALY +A (Agy 1/3 


-AALAL A, 19 -2(AZ (AZ, + (AZ 1/27 (B12) 


ALIA) FAN ERA I4) £A OTRACICA eA DIS 
HALIGZ Y HARV 14+ ALIO € (AR) 14 ALAR) + (Az) 1)/9. 


Since A? =(A’)’, the three eigenvalues A?(i=1,2,3) of A? are all real. This 


implies that 


g/4-f^/21&0. (B13) 
Let 
-g/2* dg 14 f^ 121 -|s|e" (0€ 05 7), (B14) 
we have 


-gI2-4g? IA« f? 121 2|S|e ^ . (B15) 


From Eqs. B7-9, B14 and B15, we have 


A es AL e?’ Asie ef? 4 o7 dI e 9^, (B16) 
2. 42 , ial 1/3 i03 | i2n/3 1/3 igy3 
A; =A? +e"? |S)" eo? e ur. Is enm, (B17) 
A = A? «|s[" e*^ s e, (B18) 
which are, 
A? = 4^ +2|S|'° cos((0 — 42)/3), (B19) 
A = A? «2|s|" cos((0 - 22/3), (B20) 
A? = A? «2|S[" cos(0/3). (B21) 
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Part C. Find the matrix A? from Eq. A= Üdiag (A? yu t (Eq. 14 in the main text). 


Solution: 
2 13 2 1/3 
Let AA; -2|S| cos((0—47z)/3) , AA; -2|s| cos((0—27z)/3) , and 
2 1/3 
AA; —2|S| cos(@/3), we have 
$a» 2 
A, =A ctAA, (C1) 
Dis cre: 2 
A; —A t AA, (C2) 
Dine AD. 2 
A; —AÀA + AA. (C3) 
Thus we have 
AS: ? 
A’ + AA 0 0 
Ji PE a2 2 
iag(A ) = 0 A^ c 0 . (C4) 
E) 2 
0 0 A’ +AA; 
Because 
—18cp 
C1213 515015 515€ 
T= iðcp iðcp 
U =| —5,5034 — C12S13523€ C15C53 — 5195135938 C3503 |, (C5) 
= iôcp m - iðcp 
515853 — C1281323 C155 — $,55,3053€ C13C 93 
we have 
2 2 2 
A. AL, A. 
2 2 2 2 
3 A, Apu A 
2 2 2 
Ac Au A. 
id cp 
VIE 515045 51€ 
EEN = iðcp u iðcp 
=| 75,55€»; 7 C128138238 Ci2C23 T 812813823E Ci13555 
= idcp o A iócp 
515853 — Ci5813C53€ C15854 — 8128130238 C13C 93 
2 2 C6 
A? +AA; 0 0 (C6) 
ET 2 
0 A^ * AA; 0 
"AD 2 
0 0 A‘ +AA; 
= _ id cp =o id cp 
Ci5€;5 S12C23 T €15513553€ 515353 — C1281323 
= —idcp _ = -iðcp 
5150; C15C53 T 81581355,€ C155 — $,55,3C05,€ 
iócp 
51€ C1353 C1303 


The matrix elements of A’ can be solved as 
14 


2 2 2.2 2 2 22 2 2 2 
A, = A ct co c4 AA, S044 AA; t SUAA;, (C7) 


ee 13 


Date 32 -icp 2 
Apn = C12812013C23 Z Gio Sighs ase )AA; 


2 i; 2 ió 2 (C8) 
= —l0cp —l0cp 
*(€,55,0405, T S12C 3513893 ^ JAA, CSS C AA, 
2 E -i 2 
A, = (€5555€3355 — €5€,35303€ ^ JAA, (C9) 
a EN -icp 2 “Sisk 42 
+(—Cy9819€ 3893 — $5€135:405,€ AA, + €438,3C,3€ ^ AA, 
MENU Lye iB 2 
A= pu) = (7€15515€13€53 — €15€1351355,€ ^ JAA; (C10) 
NN iScp 2 iócp A A2 
*(6,55,5014C5, —5,501451435,€ ^ DAA t €45,45, ^" AAS, 
jo A 2.2 VR 2 
AGUA (5,5€54 + €,5513553  2€,55,55,55455, COS Oop ) AA 
(C11) 
2.2 2.22 2, 22442 
HCCC + 515513554 — 2€,55,551305555, COS Opp JAA, + C55,AA;, 
jo cH 2.2 9 es 2d 2 
ue (75,5545,  €551354$54  €1531554054€ C,531551485,€ ^ DAA, 
2 2.2 2 -iô 2 ið 2, 2 2 (C12) 
EE _ —lOcp lÓcp 
* (765505435, + 515514054354 — €155,5514C5,€ + €,535,551455,€ ^ JAA + C404,5,, AA; , 
2 pABN M idee 2 
Ar, = (ALY = (655561355 — €5€13515053€ ^ DAA, 
» ^ cd idee 2 iócg A A2. 
*(7€555,5€,595, — 550,55,:05€ ^ JAAS t C,5405€ ^ AA, , (C13) 


2 ag VAN TE TA ay) 2 dócp I ibep 2 
A, OA zi 815053853 * Ci2813C23823 F C1951 98) 3C3€ C155,5814553€ )AA; 


2 2.2 2 icp 3- cops 2, 2 2 
*(7€15€53853 + 515513€53553 — Ci3515513C53€. ^ + C5515514553€ JAA; + €505455,AA; , (C14) 


PES 2 2.2.2 2 
Ax, = AX ($5555, t €45515034 — 20,55,55,4054$5, COS Opp JAA; 


TT = 12 


29 2. 23-29 2 22 2 
*(65555 + Si2813C23 t 2C 12812513023523 Cos Oc») AA; + chen AA; : (C15) 


From Eqs. C7, C11 and C15, it can be seen that the value of A2, — Abg does not 


depend on A?. 


In the case of NO (m, >m, » m,) , from Eqs. 10a-c in the main text, we have 


AA; 2—4.314x10 * eV?, (C16) 
AA; 2 —3.943x10 * eV?, (C17) 
AA; 28.257 x10 ‘eV’. (C18) 


In the case of IO (m, >m, » m,) ,by the analysis in the main text, we have 


1/3 


m? /2.- A! «2|s|" cos((6—22)/3), (C19) 


mi 12- A? «2|s|" cos(0/3), (C20) 
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1/3 


m? /2= A? «2|S|" cos(9—42)/3). 


Thus we have 
m =m? =2(A; — A>) = Als] 
= NEIK sin((z = 0) / 3), 


1/3 


[cos(0 / 3) ^ cos((2zz — 0) / 3)] 


1/3 


mi -m = 2(A; — A?) =4|s| 
= 4,/3|S|"" sin(2z —6)/3). 


[cos(0/ 3) —cos((4z — 0)/3)] 


Taking the best-fit values m;-m, =7.42x10°eV* 
m; m; -2.571x10?eV^, we have 


(m —m,) / (m — m;) —sin((z —0)/3)/sin((2z —0)/3) x 0.0289. 
From the above equation, we have 
0 =175.64°, 
and thus we have 


IS| =4.225x10“eV?. 


| 


In the case of IO, we have 


1/3 


AA? =m? /2— A? -2|s|" cos(9—2)/3) = 4.084x10 *eV?, 


1/3 


AA? = m2 12-A? -2|s|" cos(0/3) 2 4.409 x10^eV?, 


1/3 


AA? = mi 12- A? -2|S|" cos((8—47)/3) =-8.447 x10 * eV?. 


In addition, in the case of | NO, the best-fit 


(C21) 


(C22) 


(C23) 


and 


(C24) 


(C25) 


(C26) 


(C27) 


(C28) 


(C29) 


values 


6, —33.44?,0,, = 49.0°,6,, 28.57? are taken; in the case of IO, the best-fit values 


Q,, — 33.45?,0,, = 49.3?,0,, 28.61? are taken. Substituting the above values into 


the matrix elements of A?” , the results are shown in Table 1 in the main text. 
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